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Abstract 


This  report  presents  a  literature  review  on  the  modeling  of  dielectric  polymer  actuators 
(DP A).  The  principle  of  actuation  of  these  actuators  and  the  underlying  physics  of  dielectrics 
is  exposed  in  order  to  introduce  all  the  sources  of  difficulty  responsible  for  the  complex 
behaviour  of  DP  As.  A  variational  formulation  is  proposed  to  model  the  behaviour  of  DP  A, 
based  on  the  variation  of  the  dielectric  enthalpy  and  on  the  virtual  work  of  the  non¬ 
conservative  part  of  the  electric  displacement.  A  finite  element  method  is  developed,  based  on 
this  formulation.  A  reverse  mapping  requirement  is  raised  that  limits  the  type  of  element  that 
can  be  used.  An  example  of  a  possible  implementation  is  presented  for  the  case  of  a  linear 
tetrahedron.  Finally,  DP  A  modeling  challenges  are  summarized  and  recommendations  are 
made  for  future  work. 


Resume 


Ce  rapport  presente  une  revue  de  la  litterature  sur  la  modelisation  des  actionneurs  a  polymere 
dielectrique  (APD).  Le  principe  de  fonctionnement  de  ces  actionneurs,  ainsi  que  les 
phenomenes  physiques  sous-jacents  sont  exposes  afin  de  montrer  toutes  les  sources  de 
difficulte  responsables  du  comportement  complexe  de  ces  APD.  Une  formulation 
variationnelle  est  proposee  pour  modeliser  le  comportement  des  APD,  basee  sur  la  variation 
de  I’enthalpie  dielectrique  et  sur  le  travail  virtuel  de  la  partie  non  conservative  du  deplacement 
electrique.  Une  methode  element  fini  est  developpee,  basee  sur  cette  formulation.  Un 
probleme  de  reversibilite  de  la  transformation  est  leve  qui  limite  1’ application  de  cette 
methode  a  certains  elements.  Un  exemple  possible  d’ implementation  est  presente  pour  le  cas 
d’un  tetraedre  lineaire.  Enfin,  les  challenges  de  la  modelisation  d’un  APD  sont  listes  et  des 
recommandations  sont  proposees  pour  de  futurs  travaux. 
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Executive  summary 

Background 

Dielectric  Polymer  Actuators  (DP As)  that  utilize  the  Maxwell  stress  effect  consist  of  a  thin 
membrane  of  non-conductive  elastomer  sandwiched  between  elastomeric  electrodes.  When  an 
electrostatic  charge  is  applied  to  the  electrodes,  the  opposite  charges  attract,  resulting  in 
compression  of  the  elastomer  in  one  direction  and  extension  in  the  other.  Dielectric  Polymers 
present  a  non-linear  response  to  a  time  varying  voltage  due  to:  (i)  the  variable  membrane 
thickness  dependency  of  the  electric  field  and  (ii)  the  variable  area  of  the  electrodes.  Maxwell 
stress  acting  on  a  compliant  elastomer  apparently  generates  the  same  effects  as  electrostriction 
(i.  e.  the  dynamic  response  is  roughly  proportional  to  the  square  of  the  driving  voltage). 
However  it  is  clearly  a  different  mechanism.  It  was  demonstrated  in  the  literature  that  DP  As 
have  a  great  potential  as  “muscle-like  actuators”,  and  they  present  a  good  impedance  match 
with  water  so  that  they  could  make  good  low  frequency  acoustic  projectors.  However,  up  to 
now,  DP  As  have  been  built  by  trial  and  error,  using  crude  analytical  models  as  design 
guidelines.  Thus,  there  is  a  need  for  a  more  reliable  design  tool  for  developing  optimized 
DP  As.  Modeling  and  characterizing  Dielectric  Polymers  (DP)  present  a  lot  of  difficult 
challenges  since  these  materials  can  present  electric,  dielectric  and  mechanical  losses,  as  well 
non-linear  material  properties. 

Principal  Results 

This  work  presents  the  state  of  the  art  in  DPA  modeling.  It  describes  the  underlying  physics  of 
DP  As  to  explain  the  origin  of  the  major  difficulties  of  DPA  modeling  and  DP 
characterization.  A  finite  element  formulation  was  developed  which  could  be  the  starting 
point  for  developing  a  finite  element  code  for  predicting  the  behaviour  of  DP  As.  The 
background  knowledge  of  this  project  is  from  the  public  domain.  This  knowledge  has  been 
enhanced  by  the  finite  element  formulation  proposed  in  this  report. 

Significance  of  Results 

This  work  represents  a  significant  milestone  in  the  TIP  project  “Dielectric  Actuators  for 
Active/Passive  Vibration  Isolation”.  It  lays  the  foundation  for  further  modelling  work  on 
DP  As,  which  will  be  crucial  for  applications  such  as  low  frequency  active  sonar  and  active 
vibration  isolation. 

Future  Work 

Future  modelling  work  will  include  (a)  Exploring  the  effects  of  electrical  and  mechanical 
properties  on  actuation  behaviour  using  simple  models,  (b)  incorporating  a  DPA  capability  in 
the  WES  suite  of  programs  for  vibration  isolation,  and  (c)  development  of  general  finite 
element  tools  for  DP  As. 
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Sommaire 


Contexte 

Les  actionneurs  a  polymere  dielectrique  (APD),  dont  le  fonctionnement  est  base  sur  I’effet  de 
contrainte  de  Maxwell,  sont  constitues  d’une  mince  membrane  faite  d'un  elastomere  non 
conducteur  dispose  entre  des  electrodes  constituees,  elles  aussi,  d'un  elastomere.  Lorsqu’une 
charge  electrostatique  est  appliquee  sur  les  electrodes,  les  charges  opposees  s’attirent  et 
provoquent  la  compression  de  1’ elastomere  dans  une  direction  et  son  allongement  dans  1’ autre 
direction.  Les  polymeres  dielectriques  ont  une  reponse  non  lineaire  lorsqu’ils  sont  exposes  a 
une  tension  qui  varie  avec  le  temps,  et  ce,  pour  les  raisons  suivantes  :  i)  le  champ  electrique 
depend  de  I’epaisseur  de  la  membrane  et  ii)  la  surface  des  electrodes  est  variable.  La 
contrainte  de  Maxwell  qui  s'exerce  sur  un  elastomere  souple  produit  apparemment  des  effets 
identiques  a  relectrostriction  (c.-a-d.  que  la  reponse  dynamique  est  approximativement 
proportionnelle  au  carre  de  la  tension  appliquee).  Toutefois,  il  s'agit  nettement  d'un 
mecanisme  different.  Scion  la  litterature,  les  APD  presentent  un  potentiel  eleve  comme 
«  actionneurs  de  type  musculaire  »  et  pourraient  constituer  de  bons  projecteurs  acoustiques  de 
faible  frequence  en  raison  de  leur  impedance  qui  est  semblable  a  celle  de  I’eau.  Toutefois, 
jusqu’ici,  on  a  elabore  des  APD  par  approximation,  en  se  basant,  pour  la  conception,  sur  des 
modeles  analytiques  brutes.  II  y  a  done  lieu  de  disposer  d'un  outil  de  conception  plus  liable 
pour  mettre  au  point  des  APD  optimises.  La  modelisation  et  la  caracterisation  des  polymeres 
dielectriques  comportent  de  nombreux  defis  difficiles,  car  ces  materiaux  peuvent  presenter  des 
pertes  mecaniques,  dielectriques  et  electriques,  ainsi  que  des  proprietes  de  materiaux  non 
lineaires. 

Principaux  resultats 

On  presente  ici  la  fine  pointe  de  la  technologic  en  matiere  de  modelisation  des  APD.  On  decrit 
les  proprietes  physiques  sous-jacentes  des  APD  pour  expliquer  Torigine  des  principales 
difficultes  que  comportent  la  modelisation  des  ADP  et  la  caracterisation  des  polymeres 
dielectriques.  On  a  elabore  une  formulation  par  elements  finis,  qui  pourrait  constituer  le  point 
de  depart  de  la  mise  au  point  d’un  code  d'elements  finis  pour  prevoir  le  comportement  des 
APD.  Les  connaissances  necessaires  a  la  mise  en  oeuvre  de  ce  projet  sont  du  domaine  public. 
Elles  ont  ete  approfondies  grace  a  la  formulation  par  elements  finis  proposee  dans  le  present 
rapport. 

Importance  des  resultats 

Ces  travaux  constituent  une  importante  etape  du  projet  finance  dans  le  cadre  du  FIT  intitule 
« Actionneurs  dielectriques  d'isolation  active  et  passive  des  vibrations  ».  Ils  permettent 
d'etablir  une  base  a  partir  de  laquelle  seront  effectues  d'autres  travaux  de  modelisation  des 
APD,  qui  seront  d'une  importance  cruciale  pour  des  applications  comme  les  sonars  actifs  a 
basse  frequence  et  fisolation  active  des  vibrations. 
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Futurs  travaux 

Les  travaux  ulterieurs  de  modelisation  comprendront :  a)  une  etude  des  effets  des  proprietes 
electriques  et  mecaniques  sur  raction  des  actionneurs,  a  I’aide  de  modeles  simples;  b) 
rincorporation  d’une  capacite  d’APD  dans  la  suite  de  programmes  WES  sur  Fisolation  des 
vibrations;  et  c)  Felaboration  d’outils  generaux  d'elements  finis  pour  APD. 


Beslin  O.  .  2002.  Analytical  Modeling  of  Polymer  Actuators.  DRDC  Atlantic  CR  2002- 
160,  MacDonald  Dettwiler  and  Associates. 
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1.  Introduction 


The  possible  use  of  a  dielectric  polymer  (DP)  as  a  means  of  actuation  was  initially  proposed 
by  SRI  International  in  1998  [1].  The  principle  of  actuation  of  a  Dielectric  Polymer  Actuator 
(DP  A)  can  be  summarized  as  follows:  A  DP  A  is  a  condenser  filled  with  a  polymeric 
dielectric.  The  electrodes  of  the  condenser  are  mechanically  supported  by  the  dielectric 
medium  instead  of  being  maintained  by  a  rigid  frame  as  classical  condensers.  When  an 
electric  voltage  is  applied  to  the  condenser,  the  Maxwell  electrostatic  forces  acting  on  the 
electrodes  tend  to  squeeze  the  dielectric  medium.  The  electrodes  move  toward  one  another 
until  equilibrium  is  reached  between  the  elastic  stress  of  the  dielectric  medium  and  the 
Maxwell  forces.  When  the  dielectric  is  squeezed,  depending  on  its  Poisson  ratio  and 
depending  on  the  compliance  of  the  electrodes,  it  tends  to  expand  along  the  directions  of  the 
plane  of  the  condenser.  This  type  of  actuator  benefits  from  two  mechanisms  of 
“amplification”:  The  first  mechanism  is  related  to  the  compliance  of  the  dielectric  (/)  When 
the  electrodes  tend  to  merge,  the  gap  of  the  condenser  decreases  so  that  the  electric  field 
increases.  (//)  If  the  bulk  modulus  of  the  dielectric  is  high,  and  if  the  electrodes  are  compliant, 
the  dielectric  and  the  electrodes  expand,  increasing  the  surface  area  of  the  condenser.  Thus, 
both  mechanisms  tend  to  increase  the  capacity  of  the  condenser  when  a  voltage  is  applied  to 
it.  This  coupling  between  the  Maxwell  stress  and  the  elastic  stress  leads  to  a  non  linear 
mechanism  of  actuation,  yielding  a  strain  approximately  proportional  to  the  square  of  the 
driving  voltage.  In  this  description  of  the  principle  of  actuation  of  a  DPA,  the  permittivity  of 
the  dielectric  was  assumed  to  be  a  constant.  However,  in  practice,  other  effects  can  interfere 
(in  an  additive  or  subtractive  manner)  with  these  two  “amplification  mechanisms”: 

•  Electrostriction:  Le.  The  permittivity  of  the  dielectric  medium  is  strain  dependent. 

•  Quadratic  permittivity  Le.  The  permittivity  of  the  dielectric  medium  depends  on  the 
electric  field. 

These  two  non-linear  effects  add  some  complexity  to  the  development  of  a  DPA  prediction 
code,  moreover,  they  present  difficult  challenges  for  characterizing  the  mechanical  and 
dielectric  properties  of  the  dielectric  polymer. 

An  additional  source  of  difficulty  when  modeling  and/or  characterizing  dielectric  polymers  is 
losses.  Losses  can  be  classified  as  electric,  dielectric  and  mechanical. 

•  Electric  losses: 

•  Electric  losses  are  due  to  mobile  charges  in  the  dielectric,  providing  a  reminiscent 
conductivity,  they  can  also  come  from  the  electrodes  resistivity.  Moreover,  this 
resistivity  can  depend  on  the  electrodes’  expansion. 

•  Dielectric  losses: 
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•  Dielectric  losses  are  due  to  the  fact  that  the  polarization  vector  is  not  instantaneously 
following  the  electric  field  in  the  dielectric  medium.  Complex  relaxation  mechanisms 
can  occur,  as  well  as  phase  transitions. 

•  Mechanical  losses: 

•  Silicone  elastomers  can  be  considered  as  fully  elastic,  however,  acrylate  elastomers 
can  present  a  marked  viscoelastic  behavior. 

All  these  losses  contribute  to  a  temperature  increase,  which  is  another  source  of  variation  of 
the  permittivity.  In  particular,  for  polar  rubbers,  the  polarization  is  highly  temperature 
dependent,  as  rendered  by  the  Debye  model  presented  in  section  2.4.  Moreover,  thermal 
runaway  can  occur  in  actuators  made  from  elastomers  with  a  reminiscent  conductivity. 

Several  papers  [1-4]  from  SRI  International  have  experimentally  demonstrated  the  efficiency 
and  the  potential  of  DP  A  using  two  families  of  dielectric  polymers: 

•  Silicone  rubber  films  (polar  rubber  with  a  polydimethyl  siloxane  backbone) 

•  Polyacrylates  (3M  Company  VHB  acrylic  series) 

Silicone  rubbers  are  low-loss  and  soft  polymers  (0.2  to  0.4  MPa)  with  a  fast  response 
(dynamic  response  up  to  lOOOHz)  and  a  good  tear  resistance.  Polyacrylates  present  higher 
losses  and  a  higher  modulus  (0.5-0.6MPa)  with  a  slower  response  (up  to  100  Hz).  These 
papers  have  demonstrated  the  viability  of  the  concept.  However,  further  work  is  required  to 
fully  understand  the  underlying  physics  behind  these  results  in  order  to  predict  the  behavior  of 
DPA. 

A  PhD  study  [5-7]  was  conducted  at  the  Danish  Polymer  Center  focused  on  the  poly  aery  late 
3M  VHB  4910.  A  simple  analytical  model  (single  rectangular  layer  )  was  developed  to 
explain  experimental  measurements.  It  was  found  experimentally  that  the  real  part  of  the 
permittivity  drops  when  the  film  is  stretched.  They  tried  to  attribute  this  tendency  to  the 
electrostriction  effect.  However,  the  electrostrictive  behavior  of  the  VHB  4910  was  revealed 
too  weak  to  explain  properly  this  experimental  tendency.  It  was  concluded  that  the 
discrepancies  between  the  model  and  the  experiment  were  probably  due  to  the  simplistic 
model  they  used.  Particularly,  a  stretched  rectangle  is  no  longer  a  rectangle  and  some  curving 
edge  effects  have  to  be  accounted  for.  This  study  revealed  that  even  for  a  simple  experimental 
case,  the  associated  prediction  model  must  be  sufficiently  complex  to  take  into  account  the 
geometrical  effect  of  stretching. 

The  other  big  difficulty  in  DPA  modeling  is  the  availability  of  input  parameters.  As 
mentioned  previously,  the  effects  of  intrinsic  electrostriction  and  Maxwell  stress  merge  to 
generate  the  same  global  behaviour.  Thus,  it  is  difficult  to  isolate  each  effect  and  to  measure 
its  associated  coefficients,  as  the  observed  response  is  the  result  of  a  strong  coupling  between 
these  different  effects.  Guillot  et  a/.  [8]  have  proposed  a  method  to  measure  the  electrostrictive 
coefficients  of  a  polymer  film  without  being  influenced  by  Maxwell  stress  effects.  They  have 
used  this  method  to  characterize  the  electrostrictive  coefficients  of  a  polyurethane.  This 
method  should  be  tried  to  characterize  the  dielectric  of  a  DPA.  The  weakness  of  the  method. 
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however,  is  that  the  sample  to  be  tested  must  be  sandwiched  between  two  silicone  rubber 
layers.  The  stiffness  of  the  rubber  layers  is  assumed  to  be  negligible  compared  to  the  stiffness 
of  the  sample.  This  assumption  could  probably  not  be  made  when  characterizing  silicone 
rubber  films,  however,  it  could  be  possible  to  characterize  polyacrylates  with  this  method. 

In  summary,  modeling  a  DPA  and  characterizing  a  DP  present  a  lot  of  challenges.  The  aim  of 
this  report  was  to  give  a  synthesis  on  the  state  of  the  art  in  this  domain  and  to  provide  the 
basis  for  the  development  of  a  finite  element  code  for  predicting  the  behaviour  of  DPAs. 

Section  2  presents  the  underlying  mechanism  of  dielectricity  (polarization  effect).  The 
classical  Debye  model  is  presented.  The  Clausius  Mossotti  formula  is  established  and  used  to 
explain  the  electrostriction  effect  and  why  uncompressible  elastomers  present  a  low 
electrostrictive  effect. 

Section  3  presents  the  Hamilton’s  functional  of  a  DP  (and  a  pre-stressed  DP).  The  variation  of 
the  local  dielectric  enthalpy  is  calculated  to  show  the  origin  of  electrostriction  and  quadratic 
dielectricity.  Dielectric  loss  is  presented  as  the  work  of  the  non-conservative  part  of  the 
electric  displacement  for  a  variation  of  the  electric  field.  Finally,  a  finite  element  formulation 
for  modeling  Maxwell  stress  effect  in  DP  is  presented,  calculating  the  variation  of  the 
dielectric  enthalpy.  A  reverse  mapping  requirement  issue  is  raised  that  limits  the  type  of 
element  that  can  be  used.  An  example  of  possible  implementation  is  presented  for  the  case  of 
a  linear  tetrahedron.  Then,  it  is  shown  how  the  dielectric  loss  can  be  modeled  by  adding  the 
virtual  work  of  the  non-conservative  part  of  the  electric  displacement. 
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2.  Summary  of  the  theory  of  dielectrics 

In  this  section,  the  theory  of  polarization  is  presented  to  explain  the  behaviour  of  a  DP. 

2.1  Vacuum  condenser 


Figure  1:  Vacuum  condenser 


Figure  1  represents  a  vacuum  condenser  consisting  of  two  parallel  plates  (electrodes)  of 
surface  A.  The  distance  between  the  two  plates  d  is  assumed  to  be  small  compared  to  the 
dimensions  of  the  plates.  If  a  potential  V  is  imposed  on  the  upper  plate  while  the  lower  plate  is 
left  grounded,  a  surface  density  charge  +a  appears  on  the  upper  electrode  while  a  surface 
density  charge  -cr appears  on  the  lower  one.  The  charges,  give  rise  to  an  electric  fields  which 
inside  the  condenser  is  practically  homogeneous  and  directed  perpendicular  to  the  surface.  Let 
us  call  z  the  axis  perpendicular  to  the  plates  and  oriented  toward  the  up  direction.  The  electric 
field  is  aligned  along  the  z-axis  and  its  z-axis  component  is  given  by: 


The  charge  density  on  the  electrodes  a ,  and  the  electric  field  between  the  electrodes  E,  are 
classically  related  by  the  permittivity  of  the  vacuum  Sq  as  follows: 


The  total  electric  charge  maintained  at  potential  V  is  Q=oA. 
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The  capacity  of  the  condenser  C  is  such  that  Q=CV,  which  can  be  expressed  as: 

C  =  £oA 
d 

The  condenser  stores  the  total  electrostatic  energy  defined  as: 

1  1  7 

U=-QV  =  -CV^ 

'  2  2 

Increasing  the  condenser  capacity  C  allows  one  to  increase  the  storable  electrostatic  energy. 
This  can  be  achieved  by  changing  the  geometric  parameters  A  and  d  and/or  by  introducing 
between  the  electrodes  a  medium  having  a  permittivity  higher  than  the  permittivity  of  the 
vacuum.  These  kind  of  media  are  referred  to  as  dielectrics. 


2.2  Dielectric  condenser 


©  +  ©  +  ©  +  © 


V 


[-\ 

w 

w 

w 

w 

{  +  ,-}:  Free  charges 
{ @ ,  O }  •  Bound  charges 

r  i  :  Polarization 


Dielfectric  medium 


©-©-©-© 


Figure  2:  Dielectric  condenser 


Figure  2  represents  the  same  condenser  as  the  one  introduced  in  the  previous  section  but 
containing  a  dielectric  medium  between  its  electrodes  instead  of  vacuum.  If  the  same 
experiment  is  done  on  such  a  condenser,  the  electric  potential  and  electric  field  {i.e.  the 
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macroscopic  electric  field)  will  remain  the  same,  however,  the  surface  density  charge  a  will 
be  greater  than  the  one  of  the  vacuum  case.  The  increase  in  a  is  due  to  the  presence  of  electric 
dipoles  {pk}  in  the  dielectric  medium.  Depending  on  the  dielectric  medium  type,  these  dipoles 
can  originate  from  different  physical  effects  (electronic  polarization,  ionic  polarization,  polar 
molecules...).  The  average  value  of  these  electric  dipoles  over  a  small  volume  element  is 
defined  by  the  vector  polarization  P.  As  shown  in  Figure  2,  these  electric  dipoles  attract  on 
the  electrodes,  some  additional  charges,  so  called  bound  charges,  ,  (cf.  red  circled  +  and  - 
signs  in  Figure  2).  The  other  charges,  so  called  free  charges,  {cf  ordinary  +  and  -  signs  in 
Figure  2)  are  the  charges  that  would  appear  alone  on  the  electrodes  if  the  condenser  was  not 
filled  with  the  dielectric  medium.  The  sum  of  the  free  charges  and  bound  charges  is  called  the 
true  charges  and  cr  actually  represents  the  surface  density  of  these  true  charges. 

Introducing  the  permittivity  of  the  dielectric  medium  the  charge  density  cris  related  to  the 
electric  field  as  follows: 

(J 
6 

As  (j  is  higher  for  a  dielectric  medium  than  for  vacuum,  it  can  be  seen  that  the  dielectric 
permittivity  is  increased  in  the  same  proportion  as  for  a,  since  the  electric  field  remains  the 
same  in  both  vacuum  and  dielectric-filled  cases  when  the  electric  potential  is  imposed  by  an 
exterior  generator. 


2.3  Macroscopic  fields  in  dielectrics 

In  a  dielectric  medium,  it  is  convenient  to  introduce  the  electric  displacement  field  D  defined 
by: 

Equation  1  D  =  S^E  +  P 

Then,  the  basic  equations  that  govern  the  behaviour  of  electric  fields  in  dielectrics  are: 
Equation  2  E  =  —V  (j) 

Equation  3  V  x  E  =  0 

Equation  4  V  •  D  =  p 

where  (j)  is  the  electric  potential  and  p  is  the  electric  charge  density.  E  is  referred  to  as  the 
macroscopic  electric  field. 

In  the  case  of  the  dielectric  condenser  presented  in  Figure  2,  the  electric  displacement  can  be 
considered  as  constant  and  collinear  to  the  electric  field  (aligned  along  z-axis).  Then,  Equation 
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4  yields  the  simple  relation  between  the  surface  density  charge  and  the  electric  displacement: 
D=c>. 


When  the  polarization  is  collinear  to  the  electric  field  (most  cases),  Equation  1  can  be  turned 
into  a  simpler  relation: 

Equation  5  D  =  sE 

Sometimes,  the  relative  permittivity  Sy  is  used,  defined  by:  s  =  8^8^ 

Under  the  same  assumption,  the  polarization  can  be  related  to  the  electric  field  by  the  relation: 

Equation  6  P  = 

where  %  is  the  electric  susceptibility  which  is  related  to  the  permittivity  by:  8^  =  \  +  x  > 


2.4  Theoretical  models  of  the  polarization  vector 

The  polarization  vector  P  was  briefly  introduced  in  section  2.2  as  “the  average  value  of 
electric  dipoles  over  a  small  volume  element”.  This  section  will  present  a  summary  of  what 
are  the  origins  of  this  polarization  vector  in  dielectrics. 

There  is  essentially  six  ways  in  which  a  material  body  can  acquire  a  dielectric  polarization. 
These  ways  are  classified  as: 

•  Pyroelectric  effect  (built-in  permanent  polarization). 

•  Ferroelectric  effect  (built-in  reversible  polarization). 

•  Piezoelectric  effect  (stress  induced  polarization). 

•  Electronic  polarization  (induced  by  a  surrounding  electric  field  acting  on 
electrons/nucleus). 

•  Ionic  polarization  (induced  by  a  surrounding  electric  field  acting  on  ions). 

•  Dipolar  polarization  (induced  by  a  surrounding  electric  field  acting  on  polar 
molecules). 

The  material  of  interest  for  building  dielectric  polymer  actuators  must  present  a  polarization 
induced  by  an  exterior  electric  field,  which  corresponds  to  the  last  three  effects  listed  above. 
Electronic,  ionic  and  dipolar  polarizations  can  be  described  (to  the  first  order)  using  the 
polarizability  parameter  a.  If  n  is  the  number  of  induced  dipoles  per  unit  volume,  the 
polarization  is  expressed  by: 

P  =  naEi 

where  Ei  is  the  local  electric  field,  which  accounts  for  a  correction  due  to  the  additional 
electric  field  induced  by  the  polarization.  A  commonly  accepted  expression  of  the  local  field 
is: 

P 

Equation  7  E^  =  E  -{ - 

36*0 

Using  Equation  6  and  Equation  7,  the  relative  permittivity  can  be  expressed  by  the  following 
formula,  classically  referred  to  in  literature  as  the  Clausius-Mossotti  formula: 
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Equation  8 


l  +  2(^) 


3  So 


1- 


In  Equation  8,  circan  stand  for  the  electronic  polarization  {nea^,  ionic  polarization  {nta^, 
dipolar  polarization  {ndOCd)  or  a  mixture  of  them  (na  =  +  n.a.  +  ). 


Polarizabilities  are  functions  of  frequency  in  the  case  of  harmonic  applied  electric  field,  since 
the  dipoles  are  not  instantaneously  following  the  driving  electric  field.  Some  relaxation  and 
restoring  force  effects  influence  the  magnitude  and  phase  response  of  the  polarization. 
Moreover,  the  dipolar  polarization  is  greatly  dependent  on  temperature.  The  Debye-Langevin 
theory  gives  (for  low  electric  fields): 

2 

Po 

Equation  9  = - 

3kT 

where  po  is  the  electric  moment  of  the  polar  molecule,  k  is  the  Boltzmann  constant  and  T  is  the 
temperature  in  degrees  Kelvin.  Equation  9  assumes  that  the  dipoles  are  free  to  rotate  and  to 
follow  instantaneously  the  driving  electric  field.  In  practice,  additional  restoring  and  friction 
forces  must  be  considered  depending  of  the  type  of  molecules  and  geometry  (polymer 
backbone  orientation  versus  electric  field  direction). 

Due  to  the  various  origins  of  polarization  and  due  to  relaxation  processes,  it  is  evident  that 
Equation  6  (even  assuming  a  frequency  dependent  susceptibility)  is  too  simple  to  account  for 
all  the  physical  phenomena  occurring  with  time  varying  electric  field  and  dielectric  thickness. 

As  a  minimum,  a  time  relaxation  r  can  be  introduced  in  Equation  6,  accounting  for  dielectric 
losses: 

Equation  10  P  +  tP  — 


•  Digression  on  Electrostriction 

The  Clausius-Mossotti  formula  given  in  Equation  8  shows  that  the  permittivity  of  a  dielectric 
must  depend  on  the  number  of  dipoles  (and  so  molecules)  per  unit  volume.  This  means  that  a 
non-constant  volume  deformation  leads  to  changes  in  permittivity.  This  effect  of  strain 
dependent  permittivity  will  be  referred  to  as  electrostrictive  effect  in  the  following  sections. 
Moreover,  if  the  main  dielectric  property  is  due  to  dipolar  polarization,  the  permittivity 
depends  on  temperature  through  the  polarizability  ad. 

Remark:  A  formula  describing  the  electrostrictive  effect  (to  the  first  order)  can  be  derived 
from  Equation  8: 
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Equation  11 


3  ^ _ 


p 


-(6’ii  +  6*22 


+  ^33) 


where  P  = 


na 


This  last  equation  allows  us  to  understand  the  origin  of  intrinsic  electrostrictive  effect:  Any 
volume  strain  (through  the  “677+6*22+6*55”  term)  will  lead  to  a  change  in  permittivity  with  a 
sensitivity  that  will  depend  on  the  parameter  p. 
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3.  Variational  formulation 


3.1  Hamilton's  functional  of  a  dielectric  polymer 

A  dielectric  polymer  presents  dielectric  and  mechanical  losses  associated  with  internal  non¬ 
conservative  forces  (hysteresis).  However,  to  clarify,  the  equations  of  motion  of  a  dielectric 
polymer  (DP)  will  be  set  up  in  this  document  using  first  the  Hamilton  principle  of  the 
associate  conservative  case  and  then,  adding  the  virtual  work  of  the  non-conservative  forces. 

3.1 .1  Description  of  the  dielectric  medium 

Figure  3  presents  the  definition  of  the  DP  medium  domain  and  variables. 

Definitions  of  the  variables  are: 

Q  and  F  :  The  volume  and  bounding  surface,  respectively,  of  the  DP. 

(xi,  X2,  X3)  :  The  co-ordinates  of  a  DP  material  point. 

Ui :  The  displacement  in  direction  i  of  a  DP  material  point. 

(p\  The  electric  potential  in  the  DP. 

Ei :  The  electric  field  in  the  DP. 

Z)i :  The  electric  displacement  field  in  the  DP. 

Fi :  The  external  surface  force  distribution  applied  on  the  boundary  of  the  DP. 

Q  :  The  external  surface  charge  distribution  applied  on  the  boundary  of  the 
DP. 
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Figure  3.  Definition  of  the  dielectric  domain  and  fields 


3.1 .2  Hamilton's  functional 

Considering  as  the  generalized  displacement  field,  {Sij  ,  Ei }  as  the 

generalized  strain  field  and  ,  D/},  as  the  generalized  stress  field,  the  DP 
Hamilton's  functional  can  be  written  as: 

//«„  ) = I  (r  -  -  v„„ + + w„„  )jt 

h 

Where  the  energy  terms  are: 

If 

r  =  —  pii.  dEl  :The  Kinetic  energy. 


Vmech  “  “  I  •  The  mechanical  potential  energy. 
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1  (• 

^eiec  ~ - D-E-dQ. :  The  dielectric  enthalpy. 

^  Q 


Wmech  -  J  •  The  work  of  the  applied  external  surface  mechanical 

r 

loads. 


r 


3.1 .3  Variation  of  the  dielectric  enthalpy  (conservative  case) 

The  enthalpy  per  unit  volume  of  a  dielectric  material  is  given  by: 

The  variation  of  the  dielectric  enthalpy  associated  with  a  virtual  variation  of 
the  electric  field  is: 

Sg  =  --D^dE,--E,dD^ 

2  2 

If  hysteresis  phenomena  are  neglected,  the  electric  displacement  is  related  to 
the  electric  field  introducing  the  dielectric  permittivity  ^ .  In  the  general  case, 
this  permittivity  is  a  function  of  the  mechanical  strain  and  electric  fields: 

D,=^{s,,E,)E, 


The  variation  of  the  electric  displacement  is  then  given  by: 


SD,  = 


E. 

■u  J 


5E,  + 


Which  leads  to  the  following  expression  of  the  dielectric  enthalpy  variation: 


Equation  12  dg  =  -(^E^dE^  - 


y^E,j 


2  1 

E^dE^ - 

'  '  2 


E  dS: 

I  ij 


According  to  classical  literature  conventions  [8],  the  three  terms  in  the 
dielectric  enthalpy  variation  of  Equation  12  correspond  to  three  physical 
effects,  classified  respectively  by  the  following  coefficient  sets: 
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^  :  The  dielectric  permittivity  coefficient 


V®,/ 


:  The  quadratic  dielectric  coefficients 


:  The  electrostrictive  stress  coefficients 


In  the  case  of  polyacrylate  dielectric  elastomers  (3M  VHB  4910)  [6],  it  was 
shown  that  the  electrostrictive  effect  was  negligible.  More  generally,  as 
shown  previously  in  Equation  1 1,  the  change  in  permittivity  is  related  to  the 
volume  strain  in  the  polymer.  Thus,  the  more  the  polymer  is  incompressible, 
the  less  it  is  electrostrictive.  The  variation  of  the  permittivity  with  E 
(quadratic  dielectric  coefficient)  is  a  second  order  dielectric  effect. 


3.1 .4  Variation  of  the  dielectric  enthalpy  (non-conservative  case) 

For  a  mono-dimensional  problem,  considering  the  principle  of  superposition, 
the  general  linear  relation  between  the  electric  displacement  field  D{i)  and  the 
time  history  of  the  electric  field  E{i)  can  be  expressed  as[9]: 

t 

D{t)  =  +  j E{u)a{t  -  u)du 

0 

Where  a(t)  is  a  decay  function,  vanishing  with  t. 

In  this  case,  the  generalized  stress  Z)(t)  must  be  divided  into  its  conservative 
and  non-conservative  parts: 

D{t)  =  ^^E{t)  +  b{t) 

With: 


t 

D{t)  =  I  E{u)a{t  -  u)du 

0 

The  first  term,  proportional  to  E{\)  corresponds  to  the  conservative  part  of  the 
electric  displacement  so  that  it  can  be  expressed  as  the  gradient  of  a  potential 
(the  dielectric  enthalpy).  This  first  part  has  been  treated  in  the  previous 
section. 
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The  second  term  D{t)  is  related  to  the  time  history  of  the  electric  field.  This 
hysteretic  part  is  non-conservative  (the  work  of  this  stress  over  a  cycle  is  not 
zero).  Rigorously,  the  use  of  Hamilton  principle  requires  all  the  forces  to  be 
conservative.  However,  the  equation  of  motion  can  be  set  up  using  first  the 
Hamilton  principle  considering  only  the  conservative  case  and  then,  the 
virtual  work  of  the  non-conservative  forces  can  be  added  to  the  equations. 

5W  =  -dg  +  b5E 

A  numerical  implementation  of  this  non-conservative  case  is  to  consider  a 
finite  vector  {ani}{m=0  to  M)  representing  the  sampling  values  of  the  function 
a(t)  with  a  sampling  period  T. 

Then  the  relation  between  D  and  E  can  be  rewritten: 

Where  and  E^  are  the  sampled  values:  Dn=D{nT)  and  En=E(nT). 
and  where 

M 

b  =ya  E 

n  m  n-m 

m=\ 

3. 1A.  1  Measurement  of  a(t) 

The  decay  function  a(t)  is  not  easy  to  measure  using  time  domain 
experiments.  Dielectric  analyzers  provide  the  complex 
permittivity  in  frequency  domain: 


The  relations  between  a(t)  and  ^  {co)  can  be  found  classically  in 
the  literature  [9]  as: 


a{t) 


y  00 

-  j  -  L  )cos(o)t)da> 


or 


a(t) 


n 


^^"{a))s,m{cot)da) 

0 
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Those  two  last  expressions  show  that  and  are  not 
independent.  It  can  be  shown[9]  that  the  real  and  imaginary  parts 
of  the  complex  permittivity  are  interrelated  by: 


Id  -CO 


or 


^  00 

f"(®)  =  -  Jte  '(/<)  - 

n\  CO  -  n 

3.2  Hamilton's  functional  of  a  pre-stressed  dielectric 
poiymer. 

3.2.1  Definitions  of  the  pre-stressed  state  and  associated  static 
and  dynamic  fields 

If  the  DP  is  pre-stressed  by  static  electric  and/or  mechanical  charge 
distributions  F®},  the  DP  domain  and  its  bounding  surface  change.  The 
new  domain  and  boundary  surface  are  noted  Q  and  F  respectively. 

The  following  conventions  are  used: 

The  displacement  field  is  expressed  as: 

0  ,  * 
u.  =  u.  +  u. 

With: 

u\  :  The  constant  displacement  in  direction  /,  induced  by  the  static  loads 
F®}. 

Ui  :  The  additional  (time  varying)  displacement  in  direction  /,  induced  by  the 
dynamic  loads  {g  ,  F  }. 

The  strain  field  is  expressed  as: 

0  * 

Sy=Sy+Sy 

With: 

6’y®:  The  constant  strain  field,  induced  by  the  static  loads  {Q®,  F®}. 
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£//*:  The  additional  (time  varying)  strain  field,  induced  by  the  dynamic  loads 

The  electric  potential  field  is  expressed  as: 

0  * 

(p  =  (p  +(p 

with: 

(fp  :  The  constant  potential  electric  field,  induced  by  the  static  loads  F®}. 

:  The  additional  (time  varying)  potential  electric  field,  induced  by  the 
dynamic  loads  {g*,  F*}. 

The  electric  field  is  expressed  as: 

e,=e^+e: 

With; 

The  constant  electric  field,  induced  by  the  static  loads  F®}. 

E^\  The  additional  (time  varying)  electric  field,  induced  by  the  dynamic  loads 

The  electric  displacement  field  is  expressed  as: 

A=A"+A* 

With: 

Z)/^:  The  constant  electric  displacement  field,  induced  by  the  static  loads 
F®}. 

Z)/*:  The  additional  (time  varying)  electric  displacement  field,  induced  by  the 
dynamic  loads  {g*,  F*}. 

The  mechanical  stress  field  is  expressed  as: 

0  * 

(j,.  =  (j,.  +  (j,. 

u  u  u 

With: 

(7if:  The  constant  mechanical  stress  field,  induced  by  the  static  loads  {Q^, 

F®}. 
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(Ji/\  The  additional  (time  varying)  mechanical  stress  field,  induced  by  the 
dynamic  loads  {g  ,  F*}. 


3.2.2  Hamilton's  functional 

Using  the  notation  presented  in  the  previous  section,  the  Hamilton's 
functional  of  a  pre-stressed  DP  can  be  expressed  as  follow: 

^2 

//•((,,(,)  =  I  (r  -  +  (rl,  +  )dt 

h 

With; 


T  -  —  ^ p  u.u*dO. :  The  Kinetic 


energy. 


1  0  ^  0  ^ 

V^ech  “  ~  J  +  ^//)(^//  S-j)dCl :  The  DP  mechanical  potential  energy. 


Veiec  “  I (A^  +  +  E])dCl :  The  DP  dielectric  enthalpy. 


^lec  ~  I  •  The  work  of  the  applied  static  external  surface 

r* 

mechanical  loads. 


Wl^ec  ~  I  •  The  work  of  the  applied  dynamic  external  surface 

r* 

mechanical  loads. 


^eiec  ~  I (p  dT  :  The  work  of  the  applied  static  external  surface  electric 

r* 

charges. 

W^iec  =  I  Q  cp  dT  :  The  work  of  the  applied  dynamic  external  surface  electric 
r* 

charges. 
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3.3  Macroscopic  electrostriction  and  finite  element  mapping 


As  mentioned  previously,  a  dielectric  polymer  can  present  locally  a  negligible  electrostrictive 
effect  (particularly  when  its  Poisson  ratio  is  close  to  0.5).  However,  when  a  finite  element  of 
DP  is  deformed  (even  if  its  volume  is  kept  constant),  the  electric  field  inside  the  DP  element 
can  change  due  to  its  nodes’  displacements  (even  if  the  electric  potential  values  are  kept 
constant  at  these  nodes).  Hence,  the  dielectric  energy  of  a  DP  element  can  vary  with  strain, 
due  to  strain-induced  variation  of  the  electric  field.  If  is  the  enthalpy  of  a  DP  element,  and 

{ui}  its  nodal  displacements,  then  an  electrostatic  stress  tensor  must  exist  so  that: 

v.,«_ 

It  will  be  shown  in  this  section  that  the  evaluation  of  this  stress  tensor  using  a  finite  element 
method  requires  an  analytical  expression  of  the  transformation  Jacobian  matrix. 

3.3.1  Mapping 

The  problem  is  discretized  in  space  and  time: 

•  The  DP  medium  is  subdivided  into  a  finite  number  of  elements  "e". 

•  Time  is  discretized  in  samples  n  with  a  sampling  period  T. 

Figure  4  presents  a  diagram  of  the  mapping  between  the  reference  element 
(ref)  and  the  real  element  (e)  at  two  instants  t={n-\)T and  t=nT. 
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Figure  4.  Finite  element  mapping 

3.3.2  Generalized  displacement  field 

The  generalized  displacement  field  at  time  t=nT  is  expressed  on  one  element 
e  as: 


(p^^\x.,nT)  =  :  The  electric  potential  field 

(x. ,nT)  =  (x. ,  :  The  mechanical  displacement  field 

Where: 

•  is  the  value  of  the  electric  potential  field  on  the  node  k  of  the 
element  e  at  time  nT 

•  ig  the  value  of  the  mechanical  displacement  field  on  the  node  k  of 
the  element  e  at  time  nT 

•  |^  (e),A:,«  I  represents  the  set  of  the  values  of  the  mechanical  displacement 
field  at  the  nodes  {k}  of  the  element  e,  at  time  nT. 
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•  Nj^  (x. ,  }  )  is  the  interpolation  function  over  the  element  (e) 

associated  to  the  node  k,  and  expressed  as  a  function  of  the  global  co¬ 
ordinates  {Xi}.  It  must  be  noted  that  those  interpolation  functions 
expressed  in  the  global  co-ordinates  are  functions  of  the  nodal 

mechanical  displacement  set 


3.3.3  Dielectric  enthalpy 

Considering  that  the  electric  and  electric  displacement  fields  are  co-linear  and 
aligned  along  the  axis  3  of  the  global  coordinates,  the  dielectric  enthalpy  of 
the  element  e  at  time  nT  is  given  by: 

gf  =  I  -^D^^\x,,nT)E^^\x„nT)dV^'^^ 


Introducing  the  notation: 


dx-. 


The  electric  field  can  be  expressed  on  the  element  {e)  as: 
E^^\x„nT)  = 


Considering  first  the  conservative  case  {D  proportional  to  £),  the  electric 
displacement  is  given  by: 

D^^\x„nT)  =  a^E^^\x„nT) 


The  dielectric  enthalpy  can  then  be  expressed  as  a  function  of  the  generalized 
displacements 


Equation  13 


g{e)  _  ^{e);q;n  jkq  ^^^{e);k;n^^^^{e);k;n 

Where: 
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The  {  /^  }  integrals  are  symmetrical  and  can  be  calculated  on  the  reference 
element  as  presented  in  Equation  14: 


Equationl4 

j 


yiref) 


Where 


dx. 


dx. 


r{ref) 


Det{J{{x\^^’^^'  is  the  determinant  of  the  Jacobian  matrix  of  the  transformation 

{^i}^{x/}  from  the  reference  element  to  the  element  (e). 

Equation  13  can  be  rewritten  in  term  of  a  capacity  matrix  as  follows: 


g 


The  capacity  matrix  depends  on  the  position  of  the  element  nodes  at  time  step  n-1  and  on  the 
displacement  at  time  step  n: 


3.3.4  Variation  of  the  dielectric  enthalpy 

The  variation  of  the  dielectric  enthalpy  is  given  by: 

Equation  15 


-  ( {xf  } ,  {uf  *'’'}) 


where  { }  are  the  electrostatic  Maxwell  nodal  forces  given  by. 


Xf’^"’"  (  {xf*’"-' } ,  {uf })  =  -(!) 


k{e)\q\n 


sc* 

duf * 


Maxwell  electrostatic  forces  are  functions  of  the  nodal  electric  potential  as 
well  as  the  derivatives  of  the  matrix  capacity  versus  the  mechanical  nodal 
displacements 
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3.3.5  Mapping  requirements 


da'^ 


To  evaluate  — •>  the  analytical  expressions  of{ — -}  as  functions 
global  coordinates  {x/}  are  required.  Possible  finite  elements  could  be 

•  Linear  tetrahedron  (4  nodes,  (f)  for  3D  problems. 

•  Linear  Triangle  (3  nodes,  Cf)  for  2D  problems) 


of  the 
[10]: 


3. 3. 5. 1  Linear  tetrahedron 


Figure  5  presents  a  linear  tetrahedron  and  Equation  16  presents  its  associated 
interpolation  and  mapping  functions  (isoparametric  element): 
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Equation  16:  Linear  tetrahedron  interpolation  and  mapping  functions: 


A'4«,)  =  fe 


The  relation  between  the  derivatives  of  the  interpolation  functions  on  the 
reference  and  real  element  is  given  by  the  Jacobian  matrix  J: 

In  the  case  of  the  linear  tetrahedron,  Jis  given  by: 


'Jn 

-x/- 

-xf 

J  1\ 

^22 

J 22 

-x/- 

Xj'’" 

-xf 

x/’" 

-xf 

^22 

^22  _ 

-x/- 

-xy 

-xf 

where  are  the  coordinates  of  the  nodes  of  the  real  element  at  the  time 

j_  ^k:n  k:n-l  ,  k:n 

step  n:x.  =  x.  +u.  . 

The  integrals  given  in  Equation  14  can  be  evaluated  since  an  analytical 
expression  of  — -  can  be  derived  by  inverting  the  Jacobian  matrix  J. 


M"  = 


The  inverse  and  the  determinant  of  the  Jacobian  matrix  are  given  by: 


54 

Sxj 

5xj 

5xj 

54 

54 

dx^ 

5x2 

5x2 

5^2 

54 

Sxj 

5X3 

5X3 
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1 

det(J) 


•^22^^33  ~  ^^32^^23 
^31^23  ~  J  11^33 
J 21^32  ~  J 31^22 


J  13^32 

^^11^^33 

J  12^31 


J  12^33 

^^13^^31 

J 32^11 


^12^ 23 

J 2\J\3 

*^11*^22 


J  13^22 
J 23^11 
^12^21 


dct(t/)  ^11^^22^33  ^^32^^23)  ~*~ ‘^12(‘^31‘^23  ^^21^^33 )  ~*~  ^^13  (^^21^^32  ^31^22^ 


Each  Jnm  terms  are  functions  of  the  previous  positions  of  the  element  nodes  at 
time  step  n-1  :  xf and  of  the  displacement  at  current  time  step  n:  wf . 


Let  us  introduce  the  following  notations: 

F„i{xr*‘'-'}Aur"n)= 


BL  ( k™“"  uur’’})Bu  ixT'"-' } .  } ) 


dx-. 


dx-. 


Det{J{{xf*’"-^},{uf *■’"})) 


and 


The  terms  Fpp^  are  functions  of  the  real  element  nodes  positions  at  time  step  n- 
1  and  of  the  nodal  displacement  at  time  step  n. 

Integrals  /  ,  ,  are  constant  terms  as  they  are  defined  on  the  reference 


57V,. 


element.  For  the  linear  tetrahedron,  the  terms - -  are  constants  equal  to  0,1 

or-l,  as; 


dN, 


=  -1 ,  for  i=\  to  3 


dN, 


i+i 


5^. 


=  1 ,  for  i=\  to,3 


The  other  terms  are  zeros. 

Using  these  notations,  the  capacity  matrix  terms  can  be  expressed  as: 
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x; 


,{{xf {u"’"]) 


{e);k;n  ) 


da’^ 

So  that  the  terms - 77-  can  be  calculated  as: 

du’l 


dUf 


~  ^O^qkpp' 


dF, 


pp 


du: 


{{xf 


{e);k;n  ) 


3.3.6  Virtual  work  of  the  non  conservative  part  of  the  electric 
displacement 

The  non-conservative  part  of  the  electric  displacement  can  be  expressed  as  a 
weighted  summation  of  the  time  history  of  the  electric  field  from  time  step  n- 
Mto  time  step  n-\,  following  the  element  in  its  motion  with  time.  Using  a 
tetrahedral  element,  the  electric  field  is  constant  in  the  whole  element  since  it 
is  calculated  by  differentiating  linear  interpolation  functions.  So  the  non¬ 
conservative  part  of  the  electric  displacement  of  the  element  (e)  can  be 
expressed  as  a  function  of  the  element  electric  field  value  with  time: 

M 

n  m  n-m 

m=\ 


The  variation  of  the  electric  field  related  to  the  variations  of  the  nodal  electric 
potential  and  nodal  displacement  is  given  by: 

The  virtual  work  of  the  non  conservative  part  of  the  electric  displacement  on 
the  element  is  given  by: 

5W= 


This  work  can  be  divided  into  the  following  electric  and  mechanical  works; 

dW  =  dW^+SW^ 

The  electric  work  can  be  expressed  as; 
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Equation  17 


5W^  = 


where  the  integrals  are  given  by: 


Qf'"= 

y(^) 


The  mechanical  work  can  be  expressed  as: 


Equation  18 


where  the  integrals  given  by: 


Integrals  and  can  be  calculated  on  the  reference  element  using 

the  same  technique  as  the  one  presented  in  the  previous  section,  using  the 
analytical  expression  of  the  Jacobian  matrix. 

Equation  15,  Equation  17  and  Equation  18  yield  the  contribution  of  the 
electromechanical  coupling  between  the  electric  potential  field  and  the 
displacement  field  to  the  general  problem.  The  contributions  of  the  other 
terms  can  be  found  in  the  classical  finite  element  literature. 

Note:  To  simplify  meshing,  linear  tetrahedrons  can  be  assembled  to  build 
triangular  prismatic  elements,  as  shown  in  the  annex  A. 
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4.  Conclusion  and  recommendations 


This  report  has  presented  the  multiple  challenges  of  modeling  a  DPA  and  characterizing  a  DP. 
It  was  demonstrated  in  the  literature  that  DP  As  have  great  potential.  However,  up  to  now,  they 
have  been  built  by  trial  and  error,  using  crude  analytical  models  as  design  guidelines.  Some 
discrepancies  between  experimental  and  theoretically  expected  results  were  left  without  any 
clear  answer,  due  to  the  complex  behaviour  of  a  DP. 

DPA  modeling  presents  challenges  at  two  different  and  complementary  levels: 

•  The  experimental  characterization  of  DP. 

•  The  modeling  itself. 

These  challenges  are  due  to  the  complex  behaviour  of  DP  which  can  be  summarized  as 
follows: 

DPs  are  soft  materials  presenting  large  strain. 

DPs  present  a  strong  coupling  between  the  electric  field  and  the  mechanical  deformation  field. 


DPs  can  present  high  mechanical  losses  (viscoelasticity),  particularly  the  acrylic  ones. 

DPs  present  high  dielectric  losses  (hysteresis  cycle  in  D,E  diagram) 

DPs  can  present  electric  losses  (reminiscent  conductivity) 

DPs  can  present  a  small  electrostriction  effect. 

The  principle  of  actuation  of  DP  As  rely  precisely  on  behaviours  (1)  and  (2).  That  is,  they  are 
wanted  behaviours.  These  behaviours  are  a  first  source  of  complexity  for  modeling  DP  As 
since  they  require  a  non  linear  model  that  must  be  solved  in  the  time  domain. 

Behaviours  (3)  to  (5)  are  unwanted  properties  that  will  have  to  be  minimized  when  developing 
new  DP  materials.  These  behaviours  decrease  the  actuator  efficiency.  Moreover,  they  make 
difficult  the  experimental  characterization  of  DPs.  Particularly,  due  to  the  strong 
electromechanical  coupling,  it  will  be  difficult  to  separate  dielectric  and  mechanical  losses, 
making  difficult  the  measurement  of  the  dielectric  and  mechanical  relaxation  functions. 

Behaviour  (6)  can  be  a  good  property  if  its  effect  tends  to  add  the  DPA  principle  of  actuation 
(/.  e.  Maxwell  stress  +  compliance  of  the  dielectric  medium).  However  it  will  make  difficult 
the  experimental  measurement  of  the  permittivity  of  the  DP  since  the  electrostrictive  effect 
and  the  DPA  principle  of  actuation  produce  the  same  effect  (strain  proportional  to  the  square 
of  the  driving  voltage). 
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This  study  has  presented  a  possible  finite  element  method  to  model  Maxwell  stress  in  a  DPA. 
This  method  requires  the  analytical  expression  of  the  Jacobian  matrix  of  the  transformation  as 
a  function  of  the  nodes’  co-ordinates  on  the  real  elements.  It  was  shown  that  a  linear 
tetrahedral  element  could  be  used.  Even  the  dielectric  losses  could  be  calculated  (provided 
that  the  relaxation  function  of  the  DP  can  be  measured).  It  will  be  interesting  to  use  this  linear 
3D  dielectric  element  to  simulate  the  experimental  results  presented  by  Kofod  et  al.  [7].  In 
this  case,  the  predicted  strain  should  be  closer  to  the  experiment  than  the  crude  model  they 
used. 

To  build  a  reliable  finite  element  code  for  predicting  DP  As  behaviour,  it  will  be  wise  to  start 
with  a  low  loss  and  low  electrostrictive  material.  To  validate  the  proposed  element,  a  soft 
elastic  and  non-dielectric  material  could  be  used  instead  of  a  dielectric  material.  This  will 
insure  a  constant  permittivity,  however  the  strain  will  be  small  so  that  high  voltages  will  be 
required,  hopefully  a  compromise  could  be  found  between  excitation  voltage  and  material 
elasticity  so  that  a  non-linear  relation  between  excitation  and  strain  could  be  measured. 


28 


DRDC  Atlantic  CR  2002-160 


5.  References 


1.  Pelrine  Ronald  E.,  Kombluh  Roy  D.,  Joseph  Jose  P.  (1998).  Electrostriction  of 
polymer  dielectrics  with  compliant  electrodes  as  a  means  of  actuation.  Sensors  and 
Actuators  A  64,  77-85. 

2.  Kornbluh  R.,  Pelrine  R.,  Pei  Q.,  Oh  S.and  Joseph  J.  (2000).  Ultrahigh  strain 
response  of  field-actuated  elastomeric  polymers.  Smart  Structures  and  Materials 
2000:  Electroactive  Polymer  Actuators  and  Devices  (EAPAD),  Proceedings  of 
SPIE  Vol.  3987,  51-5064. 

3.  Heydt,  R.,  Pelrine  R.,  Joseph  J.,  Eckerle,  J.  and  Kombluh  R.  (2000). Acoustical 
performance  of  an  electrostrictive  polymer  film  loudspeaker.  J.  Acoust.  Soc.  Am. 
107(2),  833-839. 

4.  Pelrine  R.,  Kornbluh  R.,  Eckerle  J.,  Jeuck  P.,  Oh  S.,  Pei  Q.and  Stanford  S.  (2001). 
Dielectric  Elastomers:  Generator  Mode  Fundamentals  and  Applications.  Smart 
Stmctures  and  Materials  2000:  Electroactive  Polymer  Actuators  and  Devices 
(EAPAD),  Proceedings  of  SPIE  Vol.  4329,  148-156. 

5.  Sommer-Larsen  P.,  Hooker  J.,  Kofod  G.,  West  K.,  Benslimane  M.and  Gravesen  P. 
(2001).  Response  of  dielectric  elastomer  actuators.  Smart  Stmctures  and  Materials 
2001:  Electroactive  Polymer  Actuators  and  Devices  (EAPAD),  Proceedings  of 
SPIE  Vol.  4329,  157-156-163. 

6.  Kofod,  G.,  Kornbluh,  R.,  Pelrine,  R.  and  Sommer-Larsen,  P.  (2001).  Actuation 
response  of  poly  aery  late  dielectric  elastomers.  Smart  Stmctures  and  Materials 
2001:  Electroactive  Polymer  Actuators  and  Devices  (EAPAD),  Proceedings  of 
SPIE  Vol.  4329,  141-147. 

7.  Koffod  G.  (2001).  Dielectric  elastomer  actuators.  PhD  Thesis,  The  Technical 
University  of  Denmark. 

8.  Guillot  F.  M.and  Jarzinski  J.  (2001).  Measurement  of  electrostrictive  coefficients 
of  polymer  films.  J.  Acoust.  Soc.  Am.  1 10(6),  2980-2990. 

9.  Frohlich,  H.,(1958).  Theory  of  dielectrics  dielectric  constant  and  dielectric  loss. 

2^^  ed.  Oxford  University  Press,  London,  pi 92. 

10.  Zienkiewicz  O.  C.,  (1977).  The  finite  element  method.  3"^^  ed.  McGraw-Hill  Book 
Company  (UK)  Limited,  London,  p  787. 


DRDC  Atlantic  CR  2002-160 


29 


Annexe  A:  Example  of  tetrahedral  elements 
assembly 


Three  tetrahedral  elements  can  be  assembled  as  follows,  to  make  a  one  triangle  prismatic 
element  to  simplify  meshing  (and  post-processing  data  visualisation). 
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List  of  symbols/abbreviations/acronyms/initialisms 


DND 

DP 

DPA 


Department  of  National  Defence 
Dielectric  Polymer 
Dielectric  Polymer  Actuator 
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